Efimov Trimer Formation via Ultracold Four-body Recombination 
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We discuss the coUisional formation of Efimov trimers via ultracold four-body recombination. In 
particular, we consider the reaction A-\-A-\-A-\-B^Az-\-B with A and B ultracold atoms. We obtain 
expressions for the four-hoAj recombination rate and show that it reflects the three-hoAy Efimov 
physics either as a function of collision energy or as a function of the two-body s-wave scatter- 
ing length between A atoms. In addition, we briefly discuss issues important for experimentally 
observing this interesting and relatively unexplored process. 
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PACS numbers: 

The exquisite experimental control possible in ultra- 
cold atomic gases has enabled the observation of unique 
and bizarre quantum states. For instance, the man; 
body phenomena of Bose-Einstein condensation [l], 
and Fermi degeneracy [1, |3| have both been observed 
and used for a wide range of studies Weakly-bound 
diatomic Feshbach "halo" molecules have also been ob- 
served 0] and their dynamics probed [§]. Especially rel- 
evant to the present work, the Efimov effect Q was fi- 
nally confirmed experimentally via three-body recombi- 
nation of cesium atoms [lo| . This experimental obser- 
vation of Efimov physics was preceded by a consider- 
able amount of theoretical work due to its general im- 
portance for ultracold three-body collisions (see, for in- 
stance, Refs. [O, 12, 13, 3]) J and we now have a remark- 
ably complete characterization of these processes. 

Our goal for this Letter is to examine four-body colli- 
sions in an ultracold two-component gas. In particular, 
we want to explore the possibility of producing Efimov 
trimers via four-body recombination. Since the experi- 
mental evidence for the Efimov effect is from low-energy 
three-body collisions rather than from bound trimers, 
producing the trimers is a natural next step. One scheme 
for doing this in the tight confinement of an optical lat- 
tice has been proposed but here we investigate one 
possibility for producing them in free space. 

Besides the intrinsic interest of Efimov states, it is 
also important to develop the theory of the fundamental 
process of four-body recombination quantum mechani- 
cally. Compared with three-body collisions, our knowl- 
edge of ultracold four-body collisions is still quite rudi- 
mentary. The reason is clear: solving the Schrodinger 
equation with three additional degrees of freedom is a 
much more difficult task. There have, of course, been 
many studies of four-body systems, but only recently 
have some relevant to ultracold quantum gases begun to 
appear 3, 17, 18, l3| — none of which have addressed 



four-body recombination. 

The system we consider in this Letter is an ultracold 
mixture of atoms A and B. We take A to be bosons, 
but all we need specify now about B is that they are 
distinguishable from A. Our goal is to produce Efimov 



trimers A^ via the four-body recombination process 



A + A + A + B 



+ B. 



(1) 



To this end, we will assume that the two-body s-wave 
scattering length among A atoms gaa is infinite to give 
the most favorable case for Efimov states. We will also 
assume that the interspecies scattering length qab is fi- 
nite and that any dimer states, which are likely to be 
present in real systems, lie much deeper than the weakly 
bound Efimov trimer states. 

Our treatment of this process is asymmetrical in A and 
B, mirroring the differences in their relative scattering 
lengths. The two-body interactions for the A^ subsystem 
are well-described by the zero-range model |9|]. The ad- 
vantage of this model is that the three-body solutions are 
known and especially simple in the \aAA\ oo limit 0]. 
Note that extending this model to four-body systems 
leads to equations that are not straightforward to solve. 
Instead, we borrow an idea from Rydberg physics 13, 21 1 
— and from many-body theories of Bose-Einstein con- 
densates (BECs) as well — to model the AB inter- 
action. Since the wavelength of the A atoms, either when 
they are free or in the Efimov molecule, is much larger 
than the size of atom B, the AB interactions can be ap- 
proximated by the Fermi contact potential [23| (atomic 
units will be used throughout) 



IJ-AB 



1,2,3, 



(2) 



where ^ab is the two-body reduced mass and 7-^4 is one 
of the AB interparticle distances. The Jacobi coordi- 
nates we use to represent the internal motion of the four 
particles are shown in Fig[TJ 

Given the success of the adiabatic hyperspherical rep- 
resentation in describing the three-body continuum [111 . 
12l |. we will use it here to treat the four-body continuum 
as well. Since we want to use the known solutions for the 
A3 subsystem, we build the four-body hyperspherical co- 
ordinates from the three-body ones [2^. The four-body 
hyperradius R4 and hyperangle 04 are thus defined as 



and tanQ;4 = 



M3,4 P3 



(3) 



2 




• 3 



FIG. 1: The Jacobi coordinates for the four-body system. 
Atoms 1, 2, and 3 are identical bosons A with mass tua, 
while atom 4 is an atom of another species B with mass ms- 



Here, fi^—niA/V^, fi3^4=3mAmB/{3mA + mB), and 
/^4=^/M3M3,4- Finally, R3 denotes the three-body hyper- 
radius; and as. Delves' three-body hyperangle |24| . 

The representation we will use for the four-body prob- 
lem, however, is not fully adiabatic in that we will not 
include the interactions Vab in the adiabatic Hamilto- 
nian iJad- Rather, we will include them later as coupling 
between the channels. This choice, together with our 
definition of coordinates, permits separation of variables 
in the adiabatic equation, utilizing the known solutions 
for the A3 subsystem. Mathematically, this procedure 
begins with the four-body Schrodinger equation 

[TB^, + VAB{r 14) + VAB(r2A) + VAB{r3A) + H^] * = -B^-, 

where Tr^ is the hyperradial kinetic energy and 

H,d = Tn,+Tn,+Vi23, (4) 

which includes the hyperangular kinetic energies Tq. as 
well as all of the interactions among A atoms in V123. 
The notation U3 denotes collectively all of the three-body 
hyperangles; and O4, all remaining hyperangles for the 
four-body system. In the ultracold limit, only the zero 
total orbital angular momentum solution is relevant by 
the generalized Wigner threshold law [2^. And, since 
Efimov states only exist for zero orbital angular momen- 
tum of A3, the angular momentum of B relative to the 
trimer must also be zero. The channels for the four-body 
problem are thus defined from 



(5) 



Separation of variables allows {v = {a, n}) 

$(,'*)(i?4;^^4,r!3) = u„„(i?4;^^4)*i^)(r!3). (6) 

Because we have chosen qaa ~> 00, <i>Q (fJa) does not 
depend on R3 — otherwise this separation would only 
be approximate and the resulting equations would have 
to be solved numerically. For four-body recombination, 
we will need to find not only the A3 + B bound channels, 
but also the four-body continuum channels A+A+A+B. 



For the A3 + B channels, we use (a=0) Q 



1=1 



2 sinh(soa3 ^) 
sin(2a^'') 



2fi3Ri 



(3) 



(7) 



(8) 



where the summation is over the three possible three- 
body Jacobi sets, each with its own Delves' hyperan- 
gle a^ ^ ; so~1.0062 is a universal constant. Substituting 
Eq. (O into Eq. ©, and the result into Eq. ©, gives the 
equation for the bound trimer channels: 



9 '^0 + 4 \ rr _ 

TT-o Uon — Aq Mo„, UQn — 



X2 



2M4i?I 

The physically acceptable solution of Eq. ^ is 



(9) 



uon(i?4;a4) = A'"4cosa4sin2+*''o(a4) 
3 iso Ao„ 3 iso 



(10) 

•^On 3 2 \ 

— ; -;cos ai) 



with A^4 the normalization constant. To avoid the 
Thomas collapse [26jj, we take the simple and expedi- 
ent strategy of requiring that the three-body hyperradial 
wave functions vanish for i?3^i?o, yielding a boundary 
condition on uo„: Mon(Q!4 ^ arcsin[-\/ /i3//i4i?o/i?4])=0. 
We note that none of our conclusions will depend on 
the details of this regularization — i?o merely sets the 
scale for the features we predict. Imposing this bound- 
ary condition leads to a transcendental equation for Aon , 
and it can be shown explicitly that the Efimov trimer 
energies are recovered from /7on in the limit R^—^oo. 
We include the diagonal couphng Qj^v — ^^^^I^^ to 
get the most physical adiabatic potentials W^{R4) = 
Uu{Ri)~^^Quu{Ri)- For W^{Ri) > and R4 > i?o, 



WARi)' 



2(?i+7)-.so 



-i-f£ln 

2 TT 



/i4i?4 
/i3^0 



1 



2fi4R 



(11) 



Here, 7=0.30103 and n=l,2,3,. . . labels the trimer states. 

We follow the same logic for the four-body continuum 
channels A + A+A + B. The only difference is that 
is now a three-body continuum function, but still for 
IoaaI 00. That is, we replace sq in Eqs. ([7))- pUl) 
by isa (a>0), where Sa are real numbers determined by 



, COS(-Sa) 



8sin(-SQ). 



(12) 



The resulting Mam satisfy the same boundary condition 
at i?o as the bound channels. Asymptotically, \am 
Sa+2m+\ with m=0,l,2,. . . labeling the four-body con- 
tinuum states possible for each a. The four-body contin- 
uum potentials thus behave as I/-R4 for R4 S>i?o. 
If we expand the total four-body wave function as 



«'(i?4,f^4,f^3 



;F,(i?4)$i*'(-R4;f^4,f^3), (13) 



3 




the Landau-Zener approximation if E>W,y{R'l), where 



100 



1000 



FIG. 2: (color online) The lowest four-body adiabatic hyper- 
spherical potentials W„ multiplied by 2/i4i?| to better show 
their behavior. Black solid lines denote atom-trimer poten- 
tials for a=0 and n=l,2,3,4; blue dashed lines, four-body con- 
tinuum potentials for a=l (si =4.465) and m=0,l,2; and the 
red dotted line, the four-body continuum potential with a=2 
(s2=6.818) and m=0. We take i?o=10 a.u. for all curves. 



then the non-adiabatic couplings between four-body con- 

(3) 

tinuum and atom-trimer channels vanish since are 
independent of and form an orthonormal set. Recom- 
bination is thus driven only by the diabatic couplings 

Kv= {i^' I Vab (r 14) + Vab (r24) + Vab {r34)W) (14) 

and occurs predominantly at the crossings in Fig. [51 Fig- 
ure [5] shows the lowest atom-trimer and four-body con- 
tinuum potentials. Notice that the two sets of potentials 
cross in several places. The lowest two atom-trimer po- 
tentials, however, do not cross a continuum channel. In 
our model, then, the lowest two Efimov trimers can only 
be populated by weak non-adiabatic transitions between 
atom-trimer channels. This conclusion is independent of 
i?Oi but for more realistic, finite range two-body poten- 
tials all the channels will likely be coupled at small R4. 

Using Eq. ^ and evaluating Vi,i, numerically, we find 
that it is proportional to R^^ when R4 < |a^B|. But, 
four-body recombination should occur when R4 is com- 
parable to the size of the final three-body Efimov state. 
So, by making {uabI small compared with the size of 
the state we are interested in, we can neglect V„^. This 
requirement must be balanced against the need to have 
large off-diagonal coupling elements V^v, although there 
will be transitions so long as gab is not zero. 

In the adiabatic hyperspherical representation, four- 
body recombination Eq. ^ starts on a four-body 
continuum channel at R^—^oo. As the atoms col- 
lide, they encounter an infinite number of cross- 
ings with atom-trimer channels, one for each Efi- 
mov state, before reaching the classical turning point. 
Treating each crossing as independent from the oth- 
ers, we can estimate the transition probability using 



^4«^o./^exp 



Sn \ 



jL \c I 2j-so 



is the position 



of the crossing. In this approximation, the recombination 
probability is given by Plz = 4r(l-r) cos^Ac/). The rel- 
ative phase A(/) is approximately zero since the potentials 
are nearly parallel near the crossings, and T is 



r = exp 




27r 



KR^ (2f A£ + i) 



the 



with fc^ the wavevector and A^«y (sQ-|-2m+|)^- 
eigenvalue of Eq. ([9]) — both evaluated at the crossing. 
The unitless constant [3 originates from the evaluation 
of Vj/Vi has a weak dependence on the channel numbers, 
and is on the order of 10^^. li E < Wi,{R%), the sys- 
tem must tunnel in the initial continuum potential to 
reach the crossing and make a transition. Figure [2] and 
R% show that in this case, though, there is always an- 
other energetically accessible crossing at larger i?4 that 
will dominate the recombination. Consequently, we set 
P\,z — for all energetically closed crossings since they do 
not alter the peak structure of the total recombination 
probability. Further, we expect the first open crossing 
beyond the classical turning point to dominate all other 
open crossings since the kinetic energy grows with i?4, 
ensuring subsequent crossings will be traversed diabati- 
cally (i.e. without a transition). For the same reason, we 
neglect the possibility of transitions from an atom-trimer 
channel back to the four-body continuum at Ri>R%. 

The four-body recombination rate is related to Plz 
by /C4 c>c -^p^ where k = ^/TJZiE is the incident four-body 
wave vector. Because the couplings V^>i, are quite small, 
the peaks of Plz occur at E^W^{R'^), 



W^Rl) 



-2m- 



2Ai3,4Po 



--=^( 2n-s 



-2m-|- 



-1) 



This expression shows that for a given initial channel 
(sa, m) there is a geometrically spaced sequence of peaks 
in energy for recombination to Efimov trimers. This char- 
acteristic feature of the three-body Efimov physics thus 
comes through in the four-body physics as well. More- 
over, the spacing of the features is exp(— 27r/so) just as 
one would predict from the three-body physics. To illus- 
trate this point, we show in Fig.[3]the total recombination 
probability Pt calculated as a sum of Plz over all initial 
and final channels. Note that the width of the peaks 
increases with {qabI as does their shift from Wi,{R1). 

Up to this point, our analysis has assumed IoaaI — 00. 
If we let it be finite instead, then the above analysis ap- 
plies in the regime {uaaI ^ IoabI- To avoid the complica- 
tion of a ^2 + ^+^ continuum degenerate with the chan- 
nels of interest, we further require aAA<0 so that there 
are no weakly bound dimers. Under these conditions, the 
four-body adiabatic potentials behave as described above 
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FIG. 3: The total recombination probability Pt for atom- 
trimer channels up to n=6. (sQ,m) labels the incident chan- 
nel that contributes to the indicated peak. We have taken 
iio=10 a.u. and aAB = WO a.u.. 

in the region \aAB\ <C i?4 ^ \aAA\- For R4 > \aAA\, the 
four-body continuum potentials approach the four-body 
hyperspherical harmonic potentials Wi, ~ ^^2^1]^^^^ with 
A a non-negative integer, and the atom-trimer potentials 
approach the trimer bound energies Wi, « En . It follows 
that at energies {2ii4a\j^)~^ ■^E-^{2^4a\Q)~^ K4 keeps 
the structure described above. 

It is in the zero energy limit that the effects of finite 
o-AA reveal themselves, and we will use WKB to explore 
this limit [l^. When E ^ 0, recombination into the 
most weakly bound Efimov trimer dominates. Since the 
size of this trimer is on the order of IcaaI, we expect 
recombination to occur at i?4w|aAA|- So, following an 
analysis similar to that described in Ref. J^], we find 
that the recombination probability is 

P(x{k\aAA\ysm^{kn\aAA\+'i>), (15) 

where $ is a short-range phase independent of oaa and 
kn—^jj^exp[—{mr + 'f)/so] is the wave number for the 
final trimer state n at £^=0. When IoaaI increases by a 
factor of 22.7, a new atom-trimer channel appears and fc„ 
changes to kn+i which, in turn, changes the period of the 
ttAA-dependent oscillations. It turns out that recombina- 
tion into a particular atom-trimer state will show about 
seven full oscillations in K4 as a function of qaa- From 
the relation K4(xP/k'^, Eq. p5)) also shows that K4 will 
be constant in the threshold regime, i5<(2/X4a^^)~^. 

If the trimers cannot be experimentally observed di- 
rectly, then their production can be monitored through 
the loss of either Aor B atoms. In an ultracold mixture of 
the two, it is much better to monitor the B atoms, how- 
ever, as other few-body processes will lead to loss of A 
atoms and thus mask the effects we predict. The best sce- 
nario makes the B atoms spin-polarized fermions. In this 
case, A+A+A+B recombination is unaffected, but the 
competing loss processes involving two or more B atoms 



can be cornpletely avoided as they will be suppressed near 
threshold [25|. There remains the possibility for loss of 
B atoms in A-\-A+B collisions, but for |ayiyi|=cx) or uaa 
finite but negative, the three-body potentials are all re- 
pulsive [271] and only deeply-bound dimer channels are 
available. The various loss rates for this system are thus 
small [27!|. Monitoring the loss of B atoms should then 
provide a signature of Efimov trimer formation. 

To summarize, we have taken the first steps in under- 
standing ultracold four-body recombination into Efimov 
trimers. By carefully setting up the problem, we were 
able to obtain largely analtyical results. In the process, 
we showed that the four-body recombination should show 
prominent, geometrically spaced peaks that reflect the 
three-body Efimov physics. In fact, these peaks are sep- 
arated by precisely the factor one expects from the three- 
body physics. We also suggested one approach for exper- 
imentally determining whether four-body recombination 
was indeed occurring. Should an experiment reach this 
milestone, then our work further suggests the interesting 
possibility that one might be able to use the temperature 
to control which Efimov state is produced. 
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National Science Foundation and in part by the Air Force 
Office of Scientific Research. 
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